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I The problem of determining the electromagnetic and gravitational "self- 

o ' 

force" on a particle in a curved spacetime is investigated using an axiomatic 
approach. In the electromagnetic case, our key postulate is a "comparison 
axiom" , which states that whenever two particles of the same charge e have 
the same magnitude of acceleration, the difference in their self-force is given 
by the ordinary Lorentz force of the difference in their (suitably compared) 



X 

electromagnetic fields. We thereby derive an expression for the electromag- 
netic self-force which agrees with that of DeWitt and Brehme as corrected 
by Hobbs. Despite several important differences, our analysis of the gravita- 
tional self-force proceeds in close parallel with the electromagnetic case. In 
the gravitational case, our final expression for the (reduced order) equations 
of motion shows that the deviation from geodesic motion arises entirely from 
a "tail term" , in agreement with recent results of Mino et al. Throughout the 
paper, we take the view that "point particles" do not make sense as funda- 
mental objects, but that "point particle equations of motion" do make sense 
as means of encoding information about the motion of an extended body in 
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the limit where not only the size but also the charge and mass of the body go 
to zero at a suitable rate. Plausibility arguments for the validity of our com- 
parison axiom are given by considering the limiting behavior of the self-force 
on extended bodies. 
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I. INTRODUCTION 



In this paper, we shall investigate the motion of an isolated body coupled to classical fields 
in the limit where the spatial extent of the body is small enough that the detailed structure 
of the body is unimportant, but where the lowest order effects of the "self-field" of the body 
(which are responsible for "radiation reaction") are taken into account. Specifically, we shall 
consider (i) the motion of a charged body coupled to a Maxwell field on an arbitrary, fixed 
curved background and (ii) the motion of a massive body in an otherwise vacuum spacetime 
in general relativity. 

No difficulty of principle is encountered in the calculation of the motion of any extended 
body once one has specified what matter fields compose the body and the equations of 
motion of these matter fields. If one then gives the initial data for these matter fields (as 
well as for the classical fields to which they are coupled), the complete motion of the body 
is determined unambiguously by the full set of continuum field equations. However, in 
practice, the details of the motion of an extended body will be very complicated and will 
depend on the detailed "internal structure" of the body. Thus, if one deals with general, 
extended bodies, it is highly unlikely that any simple results can be obtained which apply 
to large classes of systems. 

One would expect the details of the internal structure of the body to become less and 
less important in the limit as one makes the body smaller. Thus, one obvious way to seek a 
class of simple, general results is to take the limit of the continuum equations as the spatial 
extent of the body goes to zero, thereby obtaining equations of motion for a "point particle" 
idealization of an extended body. However, as is well known, serious difficulties arise in 
attempting to take the point particle limit in a straightforward manner, keeping the total 
charge, e, and mass, m, of the body fixed. In the electromagnetic case, the linearity of 
Maxwell's equations allows one to make sense of the electromagnetic field of a body of finite 
charge in the point particle limit. However, the stress-energy of the electromagnetic field 
becomes singular in this limit, and the the stress-energy of the matter fields comprising the 
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body must become correspondingly singular in order for it to "hold itself together" . Hence, 
there is no reason to expect that a well defined point particle limit will exist. In general 
relativity, the situation is even worse, since the nonlinearity of Einstein's equation does not 
allow one even to make sense of the "gravitational field" of a point mass physically, 
an extended body would presumably collapse to a black hole before a point particle limit 
could be achieved. Thus, it does not appear to be mathematically or physically sensible 
to attempt to take a point particle limit — holding the charge or mass fixed — in either the 
electromagnetic or gravitational cases. 

Nevertheless, we believe that there should exist some simple and general results regarding 
the motion of bodies in the limit where the size of the body is sufficiently small (as compared, 
in particular, with the scale of variation of the background electromagnetic field and/or the 
radius of curvature of the background spacetime) and the charge and mass of the body also 
are sufficiently small that "self-field" effects do not become dominant — but are not negligible. 
Mathematically rigorous results of this sort presumably would take the form of statements 
about limits of smooth, one-parameter families of solutions to the full continuum equations 
in which both the spatial extent of the body and the total charge, e, and/or mass, m, of the 
body are simultaneously taken to zero in a suitable manner. Geodesic motion should result 
in such a limit .Q For a charged body, the Lorentz force law in the background electromagnetic 
field should then arise as a correction to geodesic motion, with the acceleration of the body 
being of order e/m. For such a charged body, we seek in this paper to find the further 
corrections to this result to order e 2 /m. Similarly, for an uncharged body, we seek to 
obtain the corrections to geodesic motion to order m in the acceleration of the body. These 
corrections arise from the "self-field" of the body and usually are referred to as "radiation 

1 This should follow directly from the theorem of Geroch and Jang Q that any world line in a 
fixed background spacetime having the property that every neighborhood of it admits a conserved 
stress energy tensor satisfying the dominant energy condition must be a timelike geodesic. 
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reaction" or "self-force" effects. In both cases, we shall not seek to obtain any additional 
corrections to the motion due to the finite size and asphericity of the bodies. There is an 
extensive literature on the equations of motion of extended bodies (see, e.g., Dixon f|] and 
references cited therein) which takes into account the lowest order deviations from Lorentz 
force or geodesic motion due to such finite size effects, but neglects the 0(e 2 /m) and 0(m) 
corrections which concern us here. Presumably, to lowest order, the combined effects of both 
types of corrections would be obtained by simply adding together the radiation reaction force 
due to self-field effects given here and "multipole" forces due to the background field which 
can be found elsewhere in the literature. 

Thus, the goal of this paper is to obtain effective equations of motion for a point particle 
which are accurate to the orders specified in the previous paragraph. It should be emphasized 
that (in contrast to some other analyses of radiation reaction phenomena) our philosophy 
is not to view a point particle of finite charge or mass as a fundamental object, but rather 
to view the point particle equations of motion we obtain as a formal device to express 
approximate results for the continuum theory, in the limit where not only the size but also 
the charge and/or mass of the body are sufficiently small. 

As already indicated above, it should be possible to derive the results we seek rigorously, 
without any additional hypotheses, by considering suitable one-parameter families of solu- 
tions in which the size, charge, and mass of the extended body go to zero in an appropriate 
manner, and suitable assumptions are made about the composition and initial state of the 
body so that its deviations from sphericity are kept under good control in this limit. How- 



ever, although we shall take a few first steps in this direction in section |II A| , we shall not 
proceed in this manner here. Rather, we shall instead attempt to "guess" the correct equa- 
tions by defining a set of axioms which we believe the total force on the particle should obey. 
These axioms uniquely determine the force, and we then shall give an explicit prescription 
which satisfies our axioms. 

Our approach bears some similarity to an approach of Penrose for obtaining the elec- 
tromagnetic self-force in the sense that, in both approaches, the self-force is given by the 



ordinary Lorentz force associated with a regularized electromagnetic field. However, the 
regularization is accomplished in very different ways in the two approaches. In the approach 
of [[| , the electromagnetic field of a point particle is regularized by means of the Kirchhoff- 
d'Adhemar formula, whereas in our approach, the electromagnetic field is effectively regu- 
larized by considering differences in the fields associated with different particle trajectories 
in (possibly different) spacetimes. We have not investigated whether these regularization 
schemes are equivalent. 

In Section [EI], we apply our approach to the case of charged point particles coupled to 
Maxwell fields on an arbitrary curved background. Our results agree with those of DeWitt 
and Brehme (as corrected by Hobbs 0), but the calculations required to obtain the 
radiation reaction force are considerably simpler and they generalize much more naturally 
to the gravitational case. In Section |ITT[ we carry out this generalization to obtain the 0(m) 
correction to geodesic motion for a point particle propagating on an otherwise vacuum 
background. Our final results agree with those recently obtained by Mino et al. |7|], but, 
again, the calculations required in our approach are considerably simpler. 

II. ELECTROMAGNETIC RADIATION REACTION 

A. Introduction and motivation for the "comparison axiom" 

In this section, we shall consider an arbitrary spacetime (M, g a b) containing an arbitrary 
timelike world line z(t) with tangent u b representing a point charge of charge e and mass m. 
We consider an arbitrary solution, F ab , of Maxwell's equations with the point particle as its 
source 



such that F is smooth except on the world line of the particle. Our aim is to give a 
prescription for assigning to each point of the world line a vector f a — which we shall refer 
to as the total electromagnetic force — such that the point particle equation of motion a b = 




(1) 



6 



f b /m, will be valid to order e 2 /m as an equation of motion for a sufficiently small, nearly 
spherical charged body, as discussed in the previous section. 

In order to motivate the axiomatic approach that we shall ultimately adopt, let us see 
what happens if we attempt to derive a formula for f a by considering a small, nearly spherical 
extended body and taking a "point particle limit" . To avoid unnecessary complications, we 
shall initially restrict attention in the discussion below to the case of a body moving in 
Minkowski spacetime. 

The first difficult issue we encounter in this program is that of obtaining a "representative 
world line" in the extended body which can be viewed as describing the motion of the body, 
so that we may contemplate the limiting behavior of this world line as the body shrinks 
toward zero spatial extent. We shall not attempt to analyze this issue here, but will merely 
assume that by methods similar to those of Beiglbock ||, it is possible to define a "center 
of mass" or "center of motion" world line such that — to an excellent approximation when 
the body is sufficiently small and nearly spherical — we have, at each point of the world line, 

p a = mu a (2) 

where u a is the tangent to the representative world line and p a is given by 

P a = I ( T body) e &c<ie, ( 3 ) 

where Tg% d is the stress-energy of the body (not including the electromagnetic field) and the 
integral is taken over the hyperplane S perpendicular to u a . Here e a bcd denotes the volume 
four- form determined by the metric, and the integrand is to be viewed as a vector- valued 
three-form; the global parallelism of flat spacetime is used here to define the integrals of 
tensor fields. It should be noted that the difficulties in controlling errors in eq. (0) and the 
corresponding equation in curved spacetime probably constitute the most serious obstacle to 
converting the arguments given in this paper into rigorous theorems about radiation reaction 
forces which do not rely on any additional axioms. 
We define the force on the extended body by 
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r = u b v b { P a ). 



(4) 



Given eq. @, we that a a = u b VbU a will be given by the projection of f a /m orthogonal to 
u a . (This projection will actually be unnecessary here in the electromagnetic case when we 
obtain the point particle equations of motion, since our final answer for f a as defined by 
eq. (|4]) will turn out to be orthogonal to u a , corresponding to m being independent of r in 
eq. (H).) We have 



where w a is the vector field which generates the map between the successive spatial 
slices £(t). Note that the value of the integral does not depend upon how we choose 
to identify successive spatial slices, i.e., it depends only on the normal component of w a 
(the lapse function) and not on its spatial projection into S(r) (the shift vector). Using the 
standard identity 



on a differential form \x (where the centered dot denotes contraction into the first index of 
the form) together with Stokes' theorem, we can rewrite the integrand to obtain 



where we now have rewritten our volume integral in more standard notation by writing 
w c dT, c in place of w c e c d e f- 

We now assume that the body is coupled to a classical electromagnetic field, F^, which 
satisfies Maxwell's equations with source given by the charge-current density, J a , of the body. 
Thus, F ab includes both the "self-field" of the body and any "background" electromagnetic 
field which may be present. We denote the stress-tensor of the Maxwell field by T^fa. For 
further generality, we also shall allow the body to be coupled to additional classical matter 




(5) 



£ w fi = w ■ djj, + d(w ■ jj) 



(6) 




(7) 



S 



(which does not couple to the electromagnetic field) with stress-energy tensor T^ t . (This 
additional matter can be thought of as a "hand" or "string" which pulls on the charged 
body. We will set T® b t = when we wish to obtain the equations of motion for a charged 
body acted upon only by an electromagnetic field, but permitting this extra coupling will 
allow us to also obtain an expression for the radiation reaction force for arbitrary world 
lines, when the charge is being "pulled".) By conservation of total stress-energy, we have 

V b [T& y + T* L + T£] = 0. (8) 

By Maxwell's equations, we have 

V b T E t = -F ab J b , (9) 
and therefore we can rewrite eq. (|7|) as 

r=f F ac j c w b di: b + f: xt (io) 



where 



/« t = - / V & T e > c d£ c (11) 

JE(t) 

For a body which accelerates with four- acceleration a b , the lapse function is given by 

N = w a u a = 1 + r 7 a 7 (12) 

where here we have extended the definition of u a over the hypersurface by parallel transport 
(i.e., u a is the unit normal to S(r)) and r 7 are the Cartesian coordinate components of the 
displacement vector on S(r) with origin at the representative world line. Thus, we obtain 

/em = f a - Text = / F«PjpN dV (13) 

where dV denotes the ordinary volume element in the Euclidean 3-space E(r) and the Greek 
indices denote components in a global inertial coordinate system. 

Let us now simplify the situation considerably further by assuming that the body is 
exactly spherically symmetric and "at rest" with respect to u a at "time" S(r), so that 
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ja _ p(j.~j u a on J](t). (However, we make no symmetry or other assumptions concerning 
the electromagnetic field.) Then we have 



/em= f r P(r)E a NdV, 

Mr) 



(14) 



where E a = F al3 up. We define q(r) to be the total charge within radius r so that 



(15) 



Then, we have 



Jem 



s(r) Airr 2 

fPq(r 



E a NdV 



S(r) 



Anr 2 



■E a N dV - 



s(t) 47rr 2 



^Dp{E a N) dV 



3(t) ^ \47rr 2 / 



q(r)E a NdV 



q(R) 



E a NdS 



g(y) 

r=i? 47ri? 2 is(r) 47rr 2 

5 3 (r)g(r)£ a AW 

<9£ 



<9r 



■AT + f%£ Q 



S(r) 



(r) 47Tr 2 



<9r 



dV, 



(16) 



where R represents the outer radius of the charged body, e = q(R) is the total charge, 
and ()n denotes the average over the sphere r = R. 

The difficulties encountered in trying to obtain a prescription for f£ M in a simple, 
straightforward manner by taking the point particle limit of the equations of motion for 
an extended body can be seen directly from eq. (|16|). Let us split the total electromagnetic 
field, F a b } into a smooth "background piece", F^\ which has a smooth limit as the size of 
the body is shrunk to zero and a piece F^ which we may view as being "due to the charge" 
itself. (Most commonly F^' would be taken to be the retarded solution with source J a .) 
When the size of the body becomes sufficiently small (in particular, much smaller than the 
scale of variation of ), the contribution of F^' to the volume integral in eq. (|T6|) becomes 
negligible, whereas the surface term straightforwardly yields the Lorentz force law 
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/£" = eF^ b u b . (17) 

since the lapse function, N, may be approximated as 1 when the body shrinks to zero size 
with E a remaining bounded. Indeed, even if the body were nonspherical and/or not perfectly 
"at rest", there would be no significant difficultly in obtaining the lowest order finite size 
corrections to (see, e.g., Dixon [[J). 

However, the situation is completely different with regard to which contributes a 
"self-force" f^M of order e 2 . If the total charge is fixed, F^' and its spatial derivatives 
become unboundedly large as the size of the body is made small, so that the integrals 
appearing in eq. fll6|) cannot easily be controlled. These integrals could still have a well 
defined limit as a result of cancellations over the different portions of the body, but the 
situation clearly is extremely delicate. In particular, small deviations from exact sphericity 
or being exactly "at rest" (as well as small corrections due to curvature when we consider 
the motion of charged bodies in curved spacetime) could easily contribute finite corrections 
to /em - Indeed, corrections to eq. (Q) itself could be large in the point particle limit. 

The singular behavior of the self-field in the point particle limit makes it very difficult 
to directly extract from (|i6|) any guidance as to what the electromagnetic force should 
be, even in the regime where one would expect this force to be essentially composition 
independent. However, inspection of eq. ( |16|) provides a possible means of dealing with 
this difficulty, and this comprises the key new idea of this paper: Suppose that, rather 
than calculating the electromagnetic force (|T6|) on a particular body, we instead attempt 
to calculate the difference in the electromagnetic force on two bodies of the same (or very 
similar) composition, which move on different world lines in (possibly different) spacetimes. 
Then, under appropriate circumstances (see below), it seems plausible that we may identify 
neighborhoods of the bodies in such a way that the difference in the electromagnetic field 
and its relevant first spatial derivatives for the two bodies can remain bounded as both of 
the bodies shrink to zero size at the same rate. If so, the volume integral contribution to 
the difference in f£ M for the two bodies (the second term on the right side of eq. fllED) will 
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go to zero. The difference in f£ M for the two bodies will then be given by a version of the 
Lorentz force law, wherein we take the difference in the electromagnetic fields, average this 
difference over the surface of the (identified) bodies as in the first term on the right side 
of eq. ( |i6|) (with N — 1), and then let the bodies shrink to zero size. Thus, if we consider 
the difference in electromagnetic forces between two bodies rather than the force on a single 
body, the "point particle limit" should be much less delicate, and, in particular, much less 
sensitive to small deviations from sphericity, etc. (provided, of course, that these deviations 
are essentially the same for both bodies). 

Indeed, even if the bodies are in different curved spacetimes, it should be possible to keep 
the difference in the electromagnetic forces on these (suitably identified) bodies under good 
control. In a curved spacetime, the above calculations would be modified in the following 
ways. First, we must use parallel transport to define the integral expression for p a . To 
do this explicitly, it is convenient to introduce an arbitrary unit vector, k a , at a point on 
the representative world line and then parallel transport k a along the worldline. For each 
r, let E(t) be the hypersurface generated by geodesies orthogonal to u a at point z{t) of 
the representative world line. For each r, we define k a on £(t) by parallel transport along 
these geodesies, thereby defining k a in a neighborhood of the entire worldline. We define 
the four-momentum p a at point z(t) by 

k a p a = Jj a (T£ dy )e bcdej (18) 

Again, we assume that a representative world line can be found for the extended body such 
that eq. ([|) holds to an excellent approximation when the body is sufficiently small and 
nearly spherical. 

If we introduce Riemann normal coordinates at point z(t) in place of the Cartesian 
coordinates of the flat spacetime derivation, the calculation which produced eq. ( |T6D from 
eq. (|5|) is essentially unchanged, with the curvature of the background introducing only a few 
corrections to eq. ([IB]). First, a new term involving Vbk a explicitly appears in the integrand 
in eq. (0) . Further corrections also result from the fact that the Riemann normal coordinate 
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components of k a are not constant. In addition, w h dT 1 \ ) deviates slightly from NdV, where 
dV is the Riemann normal coordinate volume element for the hypersurface. However, these 
corrections all decrease with the size of the particle and should become negligible in the 
point particle limit R — > 0. 

When will two bodies (possibly in different curved spacetimes) be such that the difference 
in their electromagnetic fields will be suitably bounded as their size shrinks to zero? To 
answer this question properly, we would need to carefully examine the behavior of the 
"self-field" of extended bodies as the "point particle limit" is approached. We shall not 
attempt to analyze this here. However, a good guess as to the answer to this question can 
be obtained by examining the exterior field of a point charge in a curved spacetime, and 
finding the conditions under which the difference between the fields of two such charges is — 
with a suitable identification of neighborhoods of the world lines of the particles — suitably 
bounded as one approaches a point on the world line of the particle. 

To do so, we need to study the singular behavior of the electromagnetic field of a point 
charge in curved spacetime as one approaches the world line of the point charge. As stated 
above, we are concerned only with solutions to eq. (|I|) which are singular precisely on the 
world line of the particle itself. It follows from the general theory of propagation of singular- 
ities (see theorem 26.1.1 of Hormander 0]), that all such solutions have the same singular 
behavior, i.e., the difference between any two solutions must be smooth on the world line of 
the particle itself.0 Thus, to examine the singular behavior, it suffices to focus attention on 

2 Note that this implies, in particular, that the advanced minus retarded solution is always smooth 
throughout the spacetime — including on the world line of the particle — provided only that the 
advanced and retarded solutions themselves are nonsingular off of the worldline of the particle. 
(However, as illustrated by the example of a uniformly accelerating charge in Minkowski spacetime, 
the advanced and retarded solutions need not always be nonsingular away from the worldline of 
the particle.) 
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any particular solution. When (M, g a b) is globally hyperbolic, it is convenient to examine 
the advanced and retarded solutions. The behavior of these solutions near the world line 
of the particle can be calculated by the Hadamard expansion techniques detailed by Dewitt 
and Brehme |J. The result is 

F$ v (x) = 2eg a> g\ m [r- 2 K- 1 u a 9 b + + Vv^a 2 

Z o 

- -K- 3 a a Q b ± -ktVu 6 + — K~ l u a n b R - -*T Vii! 6 c fi c 
2 3 12 6 

+ -K- 1 tt a R b c ii c + —K- 1 u a tt b R cd tt c tt d + -K~ 1 R a c b d u c tt d 
2 12 2 

- — K- 3 u a fl b R cd u c u d + -K- 3 u a R b cde u c u d fl e 

=F -K~ 2 u a R b c u c ] 
^3 J 

/±oo 
± 2V [y G a ± > „^ (Odr' + OM, (19) 

Here, u a denotes the four- velocity of the point charge at point z{r) on its world line, and x 
denotes a point sufficiently near z(t) lying on the hypersurface generated by geodesies from 
z(t) which are orthogonal to u a . The outward-directed unit tangent at z(t) to the geodesic 
passing through x is denoted by Q a , and the affine parameter of x (i.e., the distance of x from 
the world line) is denoted by r. Primed indices refer to tensors at x, while unprimed indices 
refer to tensors at z(t), and g a ' a (x, z{t)) denotes the bitensor of geodesic parallel transport. 
The quantity k is defined by k = V - M a M b V a V6cr — a a V a o", where a denotes the biscalar 
of squared geodesic distance, which plays a fundamental role in all of these expansions. 
(The normalization of a is such that a(x,z(r)) = r 2 /2.) We have set a a = u 6 VbU a and 
a a = w fe Vt,a a . The last term in the above equation is usually referred to as the "tail term" , 
and it results from the failure of Huygen's principle in curved spacetime. In that term, 
G a , a ,,{x, z(t")) denotes the advanced/retarded Green's function for the vector potential in 
the Lorentz gauge, so that 

V 6 V,Gt , - R a b Gt = -^g aa ,5(x, z), (20) 

and denotes the proper time of the point on the world line of the charged particle which 
intersects the future/past light cone of x. (The alternating sign in front of the tail term 
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integral in eq. (|T9|) merely puts the limits in the appropriate time order.) The integral in the 
tail term runs from (r ± ±e) to ±00 with the limit e — > then being taken. 3 No distributional 
component of G a , a „(x,z(r")) is encountered in the integral, and we shall assume that the 
tail term remains smooth as x approaches z(t) (as should be the case if suitable asymptotic 
conditions are placed on the world line and the spacetime). 

In addition to obvious notational differences, our expression differs from eq. (5.12) of 
Dewitt and Brehme || in three ways. First, we have added the terms that Dewitt and 
Brehme omitted due to a trivial calculational error (see Hobbs ||). Second, we have written 
the tail term in terms of the Green's function G a , a „(x, z{r")) rather than the Hadamard 
expansion term v a ' a »(x, z(t")). The latter expression gives the correct form of the tail term 
only when x is sufficiently close to z(r"), and, in general, is not even defined for large 
separations (when, in particular, there need not be a unique geodesic joining x and z{t")). 



Finally, our sign convention for the Riemann tensor is that of ||11|| , which is opposite to that 
of Dewitt and Brehme. 

In order to find the singular behavior of F^ b ,, we need to expand the coordinate com- 
ponents of g a r a and k. Both of these quantities take particularly simple forms in terms of 
Riemann normal coordinates based at z{r). We have 

9 a 'a = 5 a ' a + \r 2 WSl 5 R a ' iaS + 0(r 3 ) (21) 

and 



« = y 1 + ra a Q a + ^r 2 u a uPWn s R ail 3s + C(r 3 ). (22) 

Substituting these expansions into eq. (|19l), we have 

1 3 1 

F a'/3'( x ) = 2e[r~ 2 M [a /fi /3 , ] - -r~ 1 (a a Q a )u[ Q filpi] + -{a a VL a ) 2 u [al VL pi] + -r^a^u^ 



3 If the worldline is not complete (i.e., if it does not extend to infinite proper time in the fu- 
ture/past), then the upper limit of the tail term integral should be the maximum/minimum proper 
time values of the curve. 
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1 3 1 

1 12 1 1 

+ -U[ a '^l/3']a — -a[ a i£lpi\ ± -d[ a iupi] + — U[ a /Qp']R — -u^'R^Vl 1 
o A o IZ o 

+ -£l[ a ' -R/3'] 7 M 7 + — M[ Q ,'fi/3']-R 7 5fi 7 fi <5 + -i?[ Q /| 7 | / 3/] ( 5'U 7 fi 5 
- ^U[ a ,VLp< ] R 1 $U 1 U S + ^M[ a / J R /3 /] 75e M 7 M <5 fi e =F ^M[ Q '-R/3'] 7 M 7 ] 

± e P° V W Gi ]a „u<*"(r") dr" + 0(r) (23) 

•'Tit 



Although this formula is explicitly for the advanced/retarded solution in a globally hyper- 
bolic spacetime, as noted above, the singular behavior of F a b will be the same as in eq. fl23|) for 
any solution of Maxwell's equations with source (|l|) in a (possibly non-globally-hyperbolic) 
spacetime, provided only that F a f, is smooth away from the world line of the particle. 

From eq. (|23|), it can be seen that the divergent terms in F^j h , depend only upon the 
four- velocity and the four-acceleration of the world line at z{r). In particular, they do 
not depend upon the spacetime curvature or derivatives of the acceleration. Furthermore, 
although many of the finite terms (which do depend upon the curvature and a a ) are direction- 
dependent and thus have singular angular derivatives on the world line, the radial derivatives 
of these terms (which is all that enters the volume term in eq. (0)) are bounded. Therefore, 
it seems plausible that if we have two bodies with the same magnitude of acceleration 
at corresponding points P and P on their representative world lines and if we identify 
neighborhoods of P and P using Riemann normal coordinates, with u a aligned with u a and 
a a aligned with a a , then the singular contributions of the "self-fields" to /| M in the point 
particle limit should cancel. Thus, the difference in fg M for the two bodies in the point 
particle limit should be given by a version of the Lorentz force law wherein we average the 
difference in the electromagnetic fields over a surface of radius r as in the first term on the 
right side of eq. ([16|) (with N — 1), and then let r — > 0. This provides the motivation for 
Axiom |l] below. 
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B. The axioms 



We now are ready to state our main axiom, the motivation for which was given in the 
previous subsection. 

Electromagnetic Axiom 1 (Comparison axiom) Consider two points, P and P, each 
lying on timelike world lines in possibly different spacetimes which contain Maxwell fields F ao 
and F a b sourced by particles of charge e on the world lines. If the four- accelerations of the 
world lines at P and P have the same magnitude, and if we identify the neighborhoods of P 
and P via Riemann normal coordinates such that the four-velocities and four- accelerations 
are identifed, then the difference in the electromagnetic forces f£ M and f£ M is given by the 
limit as r — > of the Lorentz force associated with the difference of the two fields averaged 
over a sphere at geodesic distance r from the world line at P. 

fEM-f~EM = Km((F ab -F ab ))u b (24) 

Axiom [I] is a very powerful one, since it enables us to compute the difference in electro- 
magnetic force between any two particles which have the same instantaneous acceleration. 
Thus, to obtain /| M for an arbitrary trajectory in an arbitrary curved spacetime, it suf- 
fices to know /em f° r a uniformly accelerating particle — with arbitrary acceleration a a — in 
Minkowski spacetime, with the electromagnetic field chosen to be, say, the half-advanced, 
half-retarded solution. 

Let us, then, consider this special case. By symmetry, f£ M must be proportional to a a . 
If the proportionality factor were constant, such a force would correspond merely to a "mass 
renormalization" , and could be redefined away. On the other hand, such a redefinition 
would not be possible if the proportionality factor varied with acceleration. We see no 
argument from symmetry considerations alone which would forbid the presence of such 
a term. However, this spacetime, world line, and Maxwell field possess a time reversal 
symmetry about each point on the world line, which suggests that the particle always should 
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be absorbing as much electromagnetic energy as it radiates, so the electromagnetic field 
should be doing "no net work" on the particle. This, in turn, strongly suggests that /| M = 
in this case. Indeed, if we did not have f£ M = 0, the type of calculation given in section 17.2 



of Jackson |12j would show that our resulting prescription for f£ M would fail to conserve 
energy for a point particle trajectory which begins and ends in inertial motion (where, in 
this calculation, the infinite self-energy of the Coulomb field of the particle is discarded at 
the initial and final times). This motivates the following additional axiom, which agrees 
with standard claims made in textbooks (see, e.g., Jackson flT2"|): 

Electromagnetic Axiom 2 (Flat spacetime axiom) // (M, g ab ) is Minkowski space- 
time, the world line is uniformly accelerating, and F ab is the half- advanced, half-retarded 
solution, F ab = \[Fa b + F^ b ], then f a = 0at every point on the world line. 

Note that, since the advanced and retarded solutions for F ab for a uniformly accelerating 
charge in Minkowski spacetime coincide in a neighborhood of the worldline (indeed, within 
the entire "Rindler wedge" containing the worldline [13|]), it follows immediately from Ax- 
iom |l| that we also have f a = when F ab is given by the advanced solution, F^ b) or by 
the retarded solution, F~ b . Thus, we would obtain an equivalent axiom if we replaced the 
half-advanced, half-retarded solution by the advanced solution or the retarded solution. 

In the next section, we shall use Axioms [I] and ^| together with eq. ( p3[) to compute f£ M 
for an arbitrary charged particle trajectory in an arbitrary curved spacetime. 



C. The prescription 

Let P be a point on the world line of a charged particle in a curved spacetime (M, g ab ) 
containing a Maxwell field F ab satisfying eq. ([]]), where F ab is singular only on the world line 
of the particle. For simplicity we assume that (M, g ab ) is globally hyperbolic so a unique 
retarded Green's function exists; as explained at the end of this subsection, our formulas can 
easily be generalized to the non-globally-hyperbolic case. Let u a denote the four-velocity of 
the world line at P and let a a denote its acceleration at P. 
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We may view the tangent space at P as a copy of Minkowski spacetime. We shall denote 
the origin of this tangent space by P. In this Minkowski spacetime, consider a uniformly ac- 
celerating trajectory passing through P with four- velocity u a = u a and acceleration a a = a a . 
By Axiom the electromagnetic force on this uniformly accelerating Minkowski trajec- 
tory vanishes when the electromagnetic field is given by the (Minkowski) half-advanced, 
half-retarded solution. 

In (M,g a b), we write Fj£ = F ab — F~ b , where F~ b denotes the retarded solution, and it is 
assumed that F^ is smooth on the world line of the particle. Near the actual trajectory of 
the particle (in (M,g a b)), F~ b is given by eq. (^3|). On the other hand, near the uniformly 
accelerating trajectory in the tangent space at P, the (Minkowski) half-advanced, half- 
retarded solution F ab = \{F~ h + F^ h ) also is given by eq. ( p3]) except that the "tail term" 
and all of the terms involving the curvature are absent, and there is cancellation of terms 
involving a and a 2 . Axiom [l] instructs us to subtract this Minkowski retarded solution 
from F a i, (using the exponential map — or, equivalently, Riemann normal coordinates — to 
compare them), average this difference over a sphere of radius r, and then let r — > 0. The 
electromagnetic force on the particle at P is then just the Lorentz force associated with the 
resulting field. We obtain 

/em = e(F in ) ab u b + \e\a a - a 2 u a ) + \e 2 {R\u b + u a R bc u b u c ) 

+e 2 u b [ T V [b {G~) a]c 'u c ,(T')dT' (25) 

The corresponding equation of motion of a charged particle subject to no additional (i.e., 
non-electromagnetic) external forces is then simply f£ M = ma a . 

Our result ( |25| ) agrees with that of Dewitt and Brehme as corrected by Hobbs 0. 
Although we, of course, made crucial use of the Hadamard expansion for the retarded Green's 
function (p3|), no other lengthy computations were needed in our approach, since we did not 
need to compute the behavior of the electromagnetic stress-energy tensor near the world 
line of the particle. 

Note that the first term in (^) is the ordinary Lorentz force due to the incoming field. 
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The second term corresponds to the familar flat spacetime Abraham-Lorentz damping term. 
The third term is a local curvature term, whose presence is necessary to maintain conformal 
invariance of /| M . Finally, the fourth term is the so-called "tail term" resulting from the 
failure of Huygen's principle in curved spacetime. 

Due to the presence of the Abraham-Lorentz term, the equation of motion f£ M = ma a 
shares the unphysical "runaway solutions" of the ordinary flat spacetime equation of motion. 
As in the flat spacetime case, this difficulty can be resolved through the reduction of order 
technique. An exposition of the rational for this technique as well as an explanation of how 
to implement it in a general context can be found in section IV D of ]16[. To implement 
it here, we view e = e 2 /m as a "small parameter". We differentiate eq. ( p5[) (with f^ M set 
equal to ma a ) to obtain an expression for a a , and then substitute this expression back in 
eq. ([25|), neglecting terms which are higher than first order in e. We then similarly eliminate 
the terms involving a a from the right side of eq. (^5|) . The result is 

+ - — (R a b u b + u a R bc u h u c ) + —u b f T V [ \G-) a]c u c (T') dr' (26) 
3 m m J-oo 

We believe that this equation properly describes the motion of a small, nearly spherical 

charged body in a curved spacetime, taking into account the leading order effects of the 

body's "self-field". 

Inasmuch as they require the retarded solution to be singled out, expressions (p5| ) and 
( p6|) are applicable as they stand only for a particle in a globally hyperbolic spacetime. 
However, since Axiom [I] did not require global hyperbolicity, it is clear that our axioms also 
determine the electromagnetic force and equations of motion in the non-globally-hyperbolic 
case as well. Perhaps the simplest way of generalizing our formulas to the non- globally- 
hyperbolic case is as follows: If we wish to obtain f£ u at a point z(r) on the world line of a 
charged particle in a non-globally-hyperbolic spacetime, simply choose a (sufficiently small) 
globally hyperbolic neighborhood of z{r). Eqs. ( 23 ) and ( PB"D then hold at z(t), where 



and the tail term are defined in the appropriate manner, relative to that neighborhood. 
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III. GRAVITATIONAL RADIATION REACTION 



In this section, we seek to obtain the gravitational analog of our formula (p5|) above for 
the total electromagnetic force (including radiation reaction) on a charged particle, as well 
as the analog of our equation of motion ([26|) above. The latter will provide us with the lowest 
order correction to geodesic motion of a particle resulting from radiation reaction effects. In 
our approach, we shall not make any of the slow motion or post-Newtonian approximations 
common to most other treatments of gravitational radiation reaction. On the other hand, 
the applicability of our results will be limited to the motion of a small, nearly spherical 
body. 

There are many physical and mathematical similarities in the analyses of the electromag- 
netic and gravitational radiation reaction forces, and our analysis of gravitational radiation 
reaction will ultimately closely parallel that of the electromagnetic case. However, there also 
are a number of very significant differences between these two cases. We begin our analysis 
of the gravitational case by explaining in detail the nature of these differences. 

Probably the most significant difference between the electromagnetic and gravitational 
cases concerns the formulation of the question which we would like to pose. As discussed in 
detail in the Introduction, we do not view a "point particle" as a fundamental object, but, 
instead, view the "point particle limit" as a convenient mathematical means of summarizing 
results concerning the behavior of one-parameter families of extended body solutions in the 
limit where not only the size but also the charge and mass of the body go to zero in a 
suitable manner. Nevertheless, in the electromagnetic case, there is no difficulty in making 
sense of solutions to Maxwell's equations if we let the size of the body shrink to zero keeping 
its charge fixed. This enabled us to pose (and propose an answer to) the following idealized 
question in Sec |I[ Given a solution to eq. ([[]) for a Maxwell field with point particle source, 
what is the total electromagnetic force on the charged particle? On the other hand, the 
corresponding question in the gravitational case would be: Given a solution to Einstein's 
equation with a point particle source, what is the total "gravitational force" on the particle? 
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However, as already noted in the Introduction, this question makes no sense, since there is 
no notion of a solution to Einstein's equation with a "point mass" source 

A resolution of this difficulty is suggested by the fact that we are really interested in the 
case of (small) extended bodies whose self-gravity is "weak" |]Thus, it should be adequate 
to treat the gravitational effects of the body via linearized perturbation theory off of a 
background vacuum spacetime. For linear equations, there is no difficulty in making sense 
of solutions with distributional sources, so, when working with the linearized equations, it 
becomes mathematically legal to let the size of the body shrink to zero, keeping its mass 
fixed. This suggests that we pose the following question, which is directly analogous to 
the question posed in Sec ||: Let (M,g$) be a spacetime satisfying the vacuum Einstein 
equation, let z(r) be an arbitrary timelike worldline in (M, g^), and let 7^ be a solution of 
the linearized Einstein equation sourced by a particle following this worldline. What is the 
total "gravitational force" on the particle? 

Unfortunately, the above question also suffers from serious mathematical inconsistencies: 
By the linearized Bianchi identity, the linearized Einstein equation implies exact conservation 
of the stress-energy of the (linearized) source with respect to the background metric. In the 
limit where the source is a point particle, this conservation requires the world line of the 
particle to be a geodesic of the background metric. Thus, if z{r) is not chosen to be a 
geodesic of gf b , the above question makes no sense since there does not exist any solution 
7 a fe whatsoever to the linearized Einstein equation with this source. But, a knowledge of the 
total "gravitational force" only for geodesies of g^ would not be adequate for obtaining the 
self-consistent motion of the particle under the influence of its own gravitational "self-force" , 
since such a particle will deviate from geodesic motion. 

The origin of this difficulty can be understood as follows. Even for an extended body 



However, we do not wish to preclude the possibility of eventually extending our analysis to small 
bodies with strong self-gravity; see ]i~4| . 
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with very weak self-gravity, the linearized Einstein equation does not hold exactly; rather 
there are nonlinear corrections to this equation. Although these nonlinear terms make only a 
very small correction to jab, it is precisely the presence of these terms which are responsible 
for the deviations from geodesic motion. By throwing away the nonlinear terms in 7 a t,, we 
exclude from the outset the possibility that the particle fails to move a geodesic, thereby 
making it mathematically inconsistent to study departures from geodesic motion. 

To see this more explicitly, consider the exact Einstein equation for the metric g^J + 7 a ft, 
written in the form of the linearized Einstein equation for 7 a b in the Lorentz gauge, with 
the nonlinear terms in 7^ moved to the right side of the equation (in a schematic manner) 
to aid us in viewing them as an additional "source term" : 

V {0)c Vi 0) 7 afe - 2R^ c ab d % d = -16nT ab + [nonlinear terms in lab \ (27) 
V (0)a 7a b = (28) 

where % b = 7 a6 - ^g ( ° b . 

As already noted above, for a body with weak self-gravity, the matter stress-energy T ab 
should dominate the "nonlinear terms in 7 a b" . More precisely, T ab is of order m, whereas 
if there is no incoming gravitational radiation, the nonlinear terms should have magnitude 
of order m? and higher, where m denotes the mass of the body. As we shall see in more 
detail below, a knowledge of the resulting dominant 0(m) contribution to from T ab will 
suffice for determining the leading order contribution to the self-force, so we should make 
little error by dropping the nonlinear terms. However, if we do so, there are no solutions to 
eqs.(0) and fl28|) unless V^T ab = 0. 

However, a means of dealing with this difficulty is suggested by the form in which we 
have written the equations. Even when \/^T ab ^ 0, no mathematical inconsistencies occur 
in eq. ( pT|) alone when the nonlinear terms are dropped. It is only when the Lorentz gauge 
condition ( |28| ) is adjoined to this equation that inconsistencies arise. Thus, we propose to 
simply relax the Lorentz gauge condtion so that it holds only to the required accuracy, i.e., to 
0(m). (This can be ensured by simply requiring that any "incoming radiation" contributions 
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to 7 a b satisfy the Lorentz gauge condition.) The resulting system of equations should then 
have the accuracy needed to obtain the leading order contribution to the gravitational self- 
force, but should not suffer from the mathematical inconsistencies which would occur if 
the linearized Einstein equation were used to relate •y a b to T ab . We note that our viewpoint 
appears to correspond to that taken in []/]], and similar procedures for relaxing field equations 
or gauge conditions at appropriate orders also occur in many other approaches to obtaining 
self-consistent equations of motion (see, e.g., Hl5|). 

Having reformulated the equations for 7 & in this manner, we now may consider the point 
particle limit and pose the following question: Let (M,g^) be a spacetime satisfying the 
vacuum Einstein equation, let z(r) be an arbitrary timelike worldline in (M,g^), and let 
7a6 be a solution of 

V c V c7 a6 - 2R c ab d % d = -16nm J 8(x, z(r)) u\r) dr, (29) 



where 7^ = 7 a b — 7 ~ b satisfies the Lorentz gauge condition (|28|). What is the total "gravi- 
tational force" on the particle? 

Although the above question is closely analogous to the question posed at the beginning 
of Sec y, there still remain a several notable differences between the electromagnetic and 
gravitational cases. First, since we have made a linearized approximation, it is necessary 
here that 7^ be "small" compared with the background metric g^ b . No corresponding 
restriction on was necessary in the electromagnetic case. This restriction on 7^ will 
have an important bearing on the final form of the reduced order equations of motion which 
we shall obtain at the end of this section. However, it should be noted that this restriction 
on 7^ does not actually impose any physical restriction on the applicability of our results, 
since if we wished to consider a situation where the incoming, free gravitational radiation is 
"large", we could simply incorporate this radiation into the background metric gf b \ Indeed, 
there would be no (physical) loss of generality in demanding that 7^ = 0, but we choose 
not to do so, since there are a wide variety of circumstances where it is both appropriate 
and convenient to treat the incoming radiation as a linearized perturbation. 
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The second difference concerns the status of "external forces". In the electromagnetic 
case, we were free to assume that T^ t had no coupling to the electromagnetic field. However, 
in the gravitational case, it is not consistent to assume that T^ t has no gravitational coupling; 
we must include T^ t on the right side of eq. (|29|), and take into account its contributions to 
■jab- Since, ultimately, we will set T^ t = to get the equations of motion of a freely falling 
particle, this will not be relevant for our final formula for the equations of motion. However, 
in our expression for fg,, the presence of T e " b t will make a contribution to 'jab, which must be 
included. 

A third important difference concerns the gauge invariance of our results. In the electro- 
magnetic case, both the Maxwell field, F ab , and the world line, z(t), of the particle are gauge 
invariant. Most importantly, all of the information concerning the motion of the particle is 
contained in the specification of z(t). However, in the gravitational case, neither 7^ nor 
z(t) are gauge invariant, since both can be changed by diffeomorphisms. Indeed, z(t) can 
be changed arbitrarily by diffeomorphisms. Thus, the specification of z(t) alone provides 
no information about the motion of the particle. Rather, this information is encoded in the 
joint specification of both z{r) and 7^. 

Despite the above differences, our analysis of the gravitational self-force will now proceed 
in close parallel with the electromagnetic case. In order to motivate the axioms which we 
ultimately will adopt, we consider a small, nearly spherical extended body with weak self- 
gravity, so that the spacetime metric, g ab , deviates only slightly from a vacuum solution, g a ° b \ 
We seek to obtain an equation of motion for a suitable representative world line in the body, 
expressed in terms of the structures associated with the "background spacetime" (M,g a ®). 
To do so, we view the exact four-momentum density, (T^ dy )e& c de in the spacetime (M,g a b) 
from the perspective of the background spacetime (M, g^). In parallel with eq. (|IBD , at 
a point z(r) on a representative world line in the body with tangent u a , we define the 
four-momentum p a by 

k aP a = J^\T^ dy )e bcde . (30) 
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Here is the hypersurface generated by geodesies of g$ which are orthogonal (with 
respect to g^) to u a , and the vector field kffi is given the superscript "0" in order to 
emphasize that we are extending k a off of the world line by parallel transporting it with 
respect to the background metric gf b (as opposed to g a b)- As in the electromagnetic case, we 
assume that a representative world line can be chosen so that — to an excellent approximation 
when the body is sufficiently small and spherical — we have at each point of the world line 

p a = mu a (31) 

Again, the difficulties in justifying this assumption would provide one of the more formidable 
obstacles to converting the motivational arguments given here into theorems about radiation 
reaction forces. 

From the perspective of the spacetime (M,g^), the force on the body is given by 

f a = u b V { b 0) p a . (32) 

Although we would expect g a bP a p b to be constant along the world line to the order to which 
we shall work, there is no reason why gf b p a p b need be constant to this order. Equivalently, 
if we normalize u a so that g^u a u b = —1, there is no reason why the parameter m in 
eq. ( |3~TD need be constant along the curve. If g^ b p a p b f &ns to be constant, f a will fail to 
be perpendicular to u a (in the metric g^)\ we shall retain the component of f a parallel to 
u a in our formula for the gravitational force below. However, the deviation from geodesic 
motion, u 6 v£V, is determined entirely by the projection of f a perpendicular to u a in the 
metric g^ b \ i.e., we have 

mu b vf ] u a = h^\f b (33) 

where = g^f + u a Ub and all indices here are raised and lowered using gf b \ Thus, we 
ultimately will project f a perpendicular to u a when we wish to obtain the equation of motion 
of the particle. 
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The calculation of f a proceeds in parallel with the calculation in the electromagnetic 
case. Taking k a to be parallel transported (with respect to g a °J) along the world line, we 
have 

Kf a = u b vf\k aP a ) 

= J m £w{kf^ T^ dy €bcde)- (34) 

Applying the identity (|J) and using Stokes' theorem, we obtain 

Kf a = / s(o) [V 6 ^°)T b t dy + fc?> V b T£ dy ]w c e cdef ). (35) 
In the first term, we rewrite as 

V b k^ = vf ] k^ -C\ a k c , (36) 

where 

C\ a ^ \g^ d iyT lad + Vf lM - viV) (37) 

Although Vj, ^ ** will make a nonvanishing contribution to the integrand due to the back- 
ground curvature, this contribution is easily seen to vanish in our final expression for f a 
when we take the point particle limit Q, so we will drop this contribution as well as the other 
background curvature corrections mentioned in the electromagnetic derivation as they arise 
in the calculations below. 

As in the electromagnetic case, for generality, we allow the body to be coupled to addi- 
tional classical matter with stress-energy tensor T^ t . By conservation of total stress-energy, 
we have 

V b [T^ oAy + T^] = 0. (38) 



5 Note that neither V b k a nor V b k a would be negligible in the point particle limit if k a were 
defined by parallel transport with respect to g ab rather than 
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so that 

V fe T b t dy = -V b T^. (39) 
Substituting these results in eq. (p5|) , we have 

kaf a = / e(0) -k^C\ c T£ dy w d dZ d + k a f: xt (40) 

where 

k a f^ ee | s(o) -V b T^w d dE d . (41) 

We now approximate the body to be "at rest" at time E^, so that T^ dy = pu b u c , where u b 
is the unit normal (in the metric g$) to We obtain 

*«/S = K(f a - /ext) = / s(0) fciV^VV^ - uVvfVK^ (42) 

This formula corresponds to eq. ( [14] ) in the electromagnetic case, with the expression in 
parentheses playing the role of the electric field E a which appeared there. Therefore, we can 
immediately write down the gravitational analog of eq. flTEp. We obtain 

'am 



/g = m(E%N) R 



m\r) 



k(t) Airr' 



-N + r^apE^ 



dr 



dV 



+ [terms which vanish as R — > 0] (43) 



where 



E » = I w W a 7^ 7 - u^V pll a . (44) 

In this equation and in all equations henceforth, it is to be understood that all quantities 
except 7 a b refer to the background structure, and the superscript "0" will be omitted on 
the background metric and its derivative operator. As in the electromagnetic case, a "point 
particle" limit of the right side of eq. (^) cannot be taken in a straightforward manner. 
However, we can again consider the difference in f^ on two bodies of similar composition 
that move on different world lines in (possibly different) background spacetimes. In order 
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to find the conditions under which such a difference will remain bounded as R — > 0, we 
once again study the singular behavior of the exterior field of a point source in curved 
spacetime using the Hadamard expansion techniques of Dewitt and Brehme ||. Since the 
trace-reversed metric perturbation 7ab satisfies a wave equation ( p9|) very similar to the 
equation for the electromagnetic vector potential, the Hadamard expansion goes through in 
close parallel with the electromagnetic case. The covariant expansions for 7^, and V c '%/ 6 / 
(the gravitational analogs of eq. [19]) are given by 



7^(z) = 2mg a/{a g wlb) [2r 1 u a u b ±AK 2 a a u b ] 

±m G± Vallhl ,u a "{T")u h '\T") dr" + 0(r) (45) 



and 



^dlt'vi x ) = ^ m 9c' c ga'{ a g\b'\b)[-^r 2 k 1 u a u b Q c -Ar 1 k 3 a a u b u c - r x k z u a u b a c 

--K- 5 a 2 u a u b Q c + K- 5 u a u b u c a d n d + 2«T 3 aVfi c + 2K^a a a b n c 
4 a 

±-K~ 6 a 2 u a u b u c =F AK- 4 a a u b u c T ^a a a b u c =f ^a a u b a c =f ~«rV« 6 a c 
3 3 

--K- 3 R def c u a u b u d Q e u f - 2n- 3 R def a u b u c u d n e u f + 2k~ x R c ed a u h u d n e 
3 

-2K~ l R d ab n c u d u e =f 2n- 2 R c de a u b u d u e ± 2K~ 2 R d ab u c u d u e ] 

±m V c ,G± w , 6 „ W a VVV) dr" + O(r) (46) 

(In these formulas, we have normalized the advanced and retarded Green's functions G^ ba , b , 
so that they satisfy eq. (p9|) with source — 16ng aa /g bb iS(x, z).) 

Expanding g a > a and k, we find that the Riemann normal coordinate components of 
VcT^y, in the same notation as eq. (p3|), are 

V 7 '7^/ /3 /(x) = 2vn\— 2r~ 2 u a >upiVly — 4r~ 1 a (Q ,/U / 3/)U 7 / — r^ivu^ay + r _1 w a /u / 3'f2ya <5 f2 < 5 

3 3 
+6a( Q ,/Wg')U 7 'a + -u^M^/ay (a f^) — -ivw^'SXy^a fi^) 2 + Ua'Upiiiyd 

1 2 1 

+ -W a 'M/3'fi 7 '-R(5AeK'W ( 5M' : f2 A fi K — -R a i$\ t U X VL VfupiVt^i — —afUaiUpiVtryi 

2 2 
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±m I*™ V y G± a>ipil u a " (r W) dr" + 0{r) (47) 



We have verified ([|5|), (f46|) , and fl47|) using the software package MathTensor. There 
are significant discrepancies between eqs. (]45| ) and (f46l) and the corresponding equations in 



0. Specifically, (i) there is a sign discrepancy between the term ±4K~ 2 a a u b in eq. (|45|) and 
the corresponding term in eq. (2.27) of 0, (ii) many of the terms in eq. ([46]) which do not 
involve the Riemann curvature fail to appear in eq. (2.33) of 0, and (iii) there are several 
sign discrepancies between the Riemann curvature terms of eq. (f46|) and eq. (2.33) of 0. 
The sign discrepancies in the terms involving the Riemann curvature appear to arise from 



the use by [0] of inconsistent sign conventions in various formulas (the conventions of [11 
and |J both seem to appear), but we have not attempted to track down the origin of the 
other discrepancies. 

Aside from the obvious complexity introduced by the additional index structure, there 
are two important differences between the above formulas and the corresponding formulas in 
the electromagnetic case. First, as stated above, we have assumed here that the background 
spacetime is a solution of the vacuum Einstein equation, so in the gravitational case, no 
terms are present which involve the Ricci curvature. Although it would be possible to 
repeat the above analysis by perturbing off of a nonvacuum solution, the perturbations of 
the metric and background matter would become coupled at linear order, so eq. (^) no 
longer would hold, and the entire analysis would have to be redone. Second, due to the 
presence of several terms of alternating sign in the above expressions which do not depend 
upon curvature, we see that the advanced and retarded expressions for 7^ and its first 
spatial derivative for a uniformly accelerating trajectory in flat spacetime do not agree in 
a neighborhood of the worldline of the particle. In parallel with the electromagnetic case, 
the advanced and retarded solutions can be shown to be gauge equivalent within the entire 
"Rindler wedge" containing the worldline. However, unlike the electromagnetic case, the 
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analog of the Lorentz force, —Cg c ii b u c , is not gauge invariant in this case, and it differs for the 
advanced and retarded solutions even in the limit as r — > 0. Nevertheless, it can be verified 
that this difference between the forces for the advanced and retarded solutions is parallel to 
the four-velocity of the particle. Therefore, when we project the force perpendicular to the 
particle's four-velocity to produce an equation of motion, the difference will vanish and the 
situation is effectively the same as in the electromagnetic case. 

Despite the above differences, eq. ( f47|) shares the most important property of the anal- 
ogous electromagnetic expression (^), namely the divergent terms as r — > depend only 
upon the four- velocity and four-acceleration of the particle at z(t). Therefore, in direct 
analogy with Axiom p] in the electromagnetic case, we postulate the following 

Gravitational Axiom 1 (Comparison axiom) Consider two points, P and P, each ly- 
ing on timelike world lines in (possibly different) spacetimes which contain linearized metric 
perturbations sourced by particles of mass m on the world lines (see eq. ft2Q))- If the four- 
accelerations of the world lines at P and P have the same magnitude, and if we identify the 
neighborhoods of P and P via Riemann normal coordinates such that the four-velocities and 
four- accelerations are identifed, then the difference in the gravitational forces /q and f G , is 
given by the limit as r —>■ of the difference of the effective gravitational forces averaged 
over a sphere at geodesic distance r from the world line at P. 

f a G -f a G = lim (((^V a 76 C - V b7c a ) - (\v a l bc - V,7e a ))) «V (48) 

In analogy with Axiom [|, we also postulate 

Gravitational Axiom 2 (Flat spacetime axiom) // (M, g a b) is Minkowski spacetime, 
the world line is uniformly accelerating, and 7 a fc is the half-advanced, half-retarded solution, 
lab = \[lab + lab}' ^en /g = " 8 * every point on the world line. 

As noted above, since —C£ c u b u c differs for the advanced and retarded solutions, it does 
matter in this case that we use the half-advanced, half-retarded solution in this axiom, 
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rather than, say, the advanced or retarded solution, although this difference does not affect 
the projection of the force perpendicular to the worldline of the particle. 

In parallel with the electromagnetic case, the above axioms yield the following prescrip- 
tion for the gravitational force. Let (M, g ab ) be a solution of the vacuum Einstein equation 
and let 7 a b be a solution of eq. Q25D- At a point P on the particle's world line, we compare 
7 a fc with the half-advanced, half-retarded solution for a uniformly accelerating trajectory in 
the tangent space (using the exponential map to make the comparison). The gravitational 
force is then given by calculating the difference in —C bc u b u c = (|V a 7b c — Vblc a )u b u c for 
these two fields, averaging over a sphere of radius r, and letting r — > 0. If we write j ab as 
lab = l a \ + 7afe5 the resulting expression is 

fo = m^V^t - V b ( 7 in ) c >V - m 2 (yd a + ±aV) 

+m 2 u b u c f T {\v a G bca , v - V b {G-) c a albl )u a 'u v dr ! (49) 

J — oo Z 

As anticipated, this expression contains contributions to /q parallel to the four-velocity 
u a , which merely describe the effect of the metric perturbation on the normalization of u a . 
To obtain the equations of motion, we project perpendicular to u a as in Eq. (^) above. 
This yields 

« a = (^V a 7 £ " V b ( 7 in ) c a - l -u a u d V d ^ c )u b u c - ym(a» - aV) 

+mu b u c j[j l -V a G- bca , b , - V b {G-) c \, bl - Wv d G„,)// dr' (50) 

It should be noted that an Abraham-Lorentz term of the form (a a — a 2 u a ) appears in 
eq. fl50D , but with a sign opposite to that of the electromagnetic case, corresponding to the 



'antidamping" phenomenon found by Havas [17 



Finally, we apply to this equation of motion the same reduction of order techniques that 
we applied to the electromagnetic equation of motion. When we do so, the terms on the 
right side involving ma a and ma a get eliminated in favor of terms involving m^ 1 ^. However, 
our equation Q5Tj) is valid only to linear order in both m and 7^, so it is not consistent 
to keep terms involving products of these two quantities. (This contrasts strongly with the 
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electromagnetic case, we we worked only to linear order in e 2 /m, but was allowed to be as 
large as we liked, so that (e/m)F^ could be treated as being of order unity.) Consequently, 
the reduction of order procedure in this case effectively drops the Abraham-Lorentz terms, 
leaving the "tail term" as the only contribution to the "self-force"Q: 

« a = {\^ a lll ~ V b ( 7 m ) c a - \u a u d V Vy£)«V 

+mu h u c f_j\v a G^, v - V b {G-) c \, v - l -u a u d V d G- bca , b ,)u a ' u b ' dr' (51) 

This formula agrees with the results of Mino et al. 0, although they did not include 7^ in 
their expression. Not that eq. ( |51]) has a very simple interpretation: To lowest nontrivial 
order, the particle moves on a geodesic of g^} + 7^, where 7^ = 7^ + 7atf 1 and 7^ 11 = 
lab 1 ~ |flo6 7 tad i s the l as t term in equation (|45|) . 



IV. CONCLUSIONS 

In this paper, we have taken an axiomatic approach to obtain the lowest order electro- 
magnetic and gravitational "self-forces" on a small, nearly spherical body of sufficiently small 
charge and/or mass. Our final result for the total electromagnetic force on a body (possibly 
acted upon by an external electromagnetic field as well as additional, non-electromagnetic 
"external forces") is given by eq. (^). If such a body is subject only to electromagnetic 
forces, our final result for the (reduced order) equation of motion of the body is given by 
eq. (p6|). The corresponding results for the gravitational case are given by eqs. (f49[) and 

The above electromagnetic results were derived from the two axioms given in Sec. [11 B| , 
and the corresponding gravitational results were obtained from the axioms of Sec. |TJ. Al- 
though plausibility arguments in support of these axioms were given, we did not attempt 



6 This contrasts sharply with the analysis of [17] which effectively neglected the dominant "tail 
term" . 
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to prove that the electromagnetic axioms follow as a consequence of Maxwell's equations in 
curved spacetime together with conservation of total stress-energy, nor did we attempt to 
prove that the gravitational axioms follow from Einstein's equation. Nevertheless, we believe 
that our plausibility arguments have provided the first steps in that direction. In any case, 
the problem of providing a rigorous justification for the electromagnetic and gravitational 
self-forces and the corresponding equations of motion has been reduced to the problem of 
providing a rigorous justification of our axioms. 

Finally, we note that on account of the "tail term" , our equations of motion in both the 
electromagnetic and gravitational cases are integro-differential equations which, in principle, 
require us to know the entire past history of the particle. However, if the curvature of 
spacetime is sufficiently small and the motion of the body is sufficiently "slow", one would 
expect the "tail term" to become effectively local, since the contributions to the "tail term" 
arising from portions of the orbit distant from the present position of the particle should 
become negligible. Indeed, if the tail term becomes effectively local, arguments using the 
standard dipole formula for radiated energy together with conservation of energy (see W2 



and the analogous gravitational calculation given below) suggest that if F$ = and R a b = 0, 
the "tail term" of eq. (25) should reduce to the familiar Abraham-Lorentz damping force 



; ™ _ 3 e ^- (52) 



That this is indeed the case was established by DeWitt and DeWitt [T8|] for a charged particle 
in a slow, circular orbit in linearized Schwarzschild spacetime. It should bve emphasized that 
the "true" Abraham-Lorentz force actually vanishes here — since a a = a a = for geodesic 
motion — but, remarkably, the tail term mocks up an effective Abraham-Lorentz term. 

In the gravitational case, the standard "quadrupole formula" for radiated power in the 
slow motion, weak field limit is 

where i and j are spatial indices = 1,2,3) for a set of global inertial coordinates and 
is the traceless quadrupole moment 
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Qij = q ij _ l q5 v (54) 

qij = 3 J T oo x i x j d 3 x _ ( 55 ) 

so, for a point particle, 

Q ij = 3m (W - -a^x*.^'") . (56) 



3 

By conservation of energy, we should have 

1 rd 3 Q l U 3 Q 



dt 



1 pcFQVdQij 



45 J dt 5 dt 
2 r ( d 5 Q ij 



-m 



Xj V{ dt. (57) 



15 J V dt 5 

which suggests that the radiation reaction force should be given by (see section 36.8 of 
Misner, Thorne, and Wheeler ]TT|) 



JG 15 dt 5 3 y ' 

It would be interesting to perform the gravitational analog of the analysis of DeWitt and 



DeWitt [18] to see if this formula does, indeed, arise from the tail term of eq. (pTD in the 



slow motion, weak field limit. 
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